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1 Introduction 



A project to numerically study spatial correlators G(r) of the scalar and pseudoscalar oper- 
ators Tr F 2 and TV FF in finite temperature gauge field theory has been initiated in [TJ [2] . 
Analytic computations of this correlator in next-to-leading order QCD perturbation theory 
have been carried out in [31 0]. The purpose of this note is to compute this imaginary-time inte- 
grated entirely static correlator in supersymmetric N = 4 Yang-Mills theory using AdS/CFT 
duality. There is a standard framework for this [51 El El |9] , but the computation of full 
r dependence involves some subtleties in the subtraction of divergences in the u, k plane so 
that it is perhaps motivated to report on the details of the computation. 
Concretely, we want to compute the finite T correlator, of dimension 7, 



G(r)= f dr(F\r, X )F\0,O)) T = f drG(r,x), F 2 = \F*,F*,, 
Jo Jo 



|x|- (1) 



This correlator is purely static, ui = 0, and we shall determine it by first evaluating its 3d 
Fourier transform G(k), of dimension 4, and then transforming to coordinate space: 

G( r ) = -^t dkksin(rk)G(k). (2) 



2vr 2 r jo 

Hereby it is essential to separate the vacuum part from the finite T part: 

G(k) = G vac (k) + [G(k) - G vac (k)] = G vac (k) + G T {k). (3) 

We shall include in G vac (k) all constant terms and evaluate them analytically. Transforming 
to coordinate space they lead to contact terms proportional to <5(x), but they have to be 
subtracted in numerics and their precise value is important. Furthermore, we show that at 
large k, Gr(k) can be expanded in powers of 1/k 4 and evaluate analytically the terms up to 
order 1/k 8 . 

Although the main result can be obtained by taking u = 0, it is also of interest to see 
with what happens when the real time frequency to is included. Computing G(cj,k) and, in 
particular, the spectral function p(w,k) = ImGi?(w,k) using the methods of [51 El El [8] , the 
coordinate space Green's function is given by 



/ d 3 k ix . k f°° du . , . cosha/3-rju; 

Gr ' X = / 7^3 e / — ■ I la ■ 4 

J (27:)* J it sinhj^w 
Integrating this over r we see that the static correlator is also obtained by evaluating 

GW i r dkksWk ) rt±^L. (5) 



2vr 2 ? 



7T 0J 



We shall show that the same result is obtained via ^ and ([5]), though the path via ^ is much 
longer. An essential part of the computation is again correct identification and elimination of 
divergences. Furthermore, interesting special nonstatic cases are the zero momentum temporal 
correlator G(t, k = 0) and the equal time correlator G(t = 0,r). We shall comment on these 
at the end. An AdS/CFT computation of G(t, r) has been presented in [2]. 
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2 Equations 

In the standard framework (5j El [7J |8] and notation, one takes the background 

ds 2 = b 2 (z)(-f(z)dt 2 + dx. 2 + ^) , (6) 



/(*) 

C z 4 C? N 2 

and the scalar field equation therein: 

•• (3b A ■ /w 2 fc 2 \ 

where = (f>(z,u,k) = <f>(z,K), k = |k| . Scaling all dimensionful quantities with 27, the 
equation becomes 



3 + z 



1 



r 2(1 - z 4 

where 



+ 



w 2 k 2 



(1-z 4 ) 2 1-z 4 



0, (9) 



to k 
io^z h u=—, k^z h k=—. (10) 



Since 



the Wronskian of two any linearly independent solutions of (|8j) is, integrating W/W = —P, 

where Wo is a z independent constant (but will depend on to, k). 
2.1 The static case, u = 

As explained in the introduction, our goal is to compute the correlator in the static limit 
to = 0. In this case it is a purely Euclidean quantity, and the equation to be solved is simply 

3 + z 4 ■ k 2 

r<P = 0, (13) 



Y z{\-z 4 y 1-z 4 

with the boundary conditions 0(0) = 1 and |0(1)| < 00. While the full solution of this equation 
is very complicated, its behavior at large k can be extracted. The leading term will give the 
vacuum part of the correlator, which diverges as k 4 and has to be subtracted before Fourier 
transforming our numerical results into coordinate space. The subleading terms will enable 
us to analytically compare the short-distance limit of our result with that in perturbative 
QCD [4J, and also to have better control over the numerics. 
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To find the k — >• oo limit of (13), the simplest method is to scale the variable as y = kz 
and expand in 1/k 4 . The resulting differential equations are then solved order by order in 
k~ 4 , with the requirement that the solution stay finite at z ~ 1. The leading terms and their 
expansion at small z are 

<t>(z, k) = \ y 2 K 2 (y) + £ • i [y e K 2 (y) - y 7 K^y))] + 0(k^) 

w 1 - i (A;^) 2 - i A: 4 z 4 (ln ifc + lnz + lE - \ ) + i z 4 + 0(z 6 ). (14) 

The constant term 1/10 will be important in numerics. 

This method, while intuitive, is hard to implement beyond this point, as each new order 
requires solving an inhomogeneous differential equation on the whole interval < y < oo to 
keep track of the boundary conditions. In practice it is better to follow the method of Olver 
[SI [TUJ. First we take u = z 2 as a new variable and remove the first derivative by writing 
4>{u) = W(u)yfu/(1 — u 2 ). The resulting equation for W(u) reads 



W"(u) 



+ 



3 - 6u 2 - u 4 



n(l - u 2 ) Au 2 (l-u 2 ) 2 
where, as always when the variable u is used, k = k/(2irT). 



W(u), 



(15) 



Following |10j . we then rewrite Eq. (15) in terms of £ and itf(C) given by 



c 



1/2 



dt 



and obtain the equation 



-1/2 



w"(C) 



e _3_ 

4C + 4C 2+ C 



MC), 



where 



^(C) 1280 £ 3 665600 ^ + 20" I Son 



.I,,3 
20 " 



73 u 5 + 0(« 7 ). 



The solution for w(C) can then be written as a series in inverse powers of k 2 , 

oo 

;(C) = C 1/2 ^(K 1/2 )£ 



w( 



s=0 



w(C) (16) 

(17) 
(18) 

(19) 



s=0 



where Aq = 2k 2 and the other functions A s and B s are found using the recursion relations 



1 /" c dv r 



r/)(u)A s (u) 



B a (C) = -4(0 + ^ 

A s+ i(C) = 2B S (() - (B' s (() + f dv^(v)B s (v). 

Jo 



(20) 
(21) 



The leading term, given by Aq alone, reproduces eq. (14). For the following terms we note 
that for the purposes of computing the correlator, it is sufficient to work out the power series 
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of 4>{u) around u = to order u 2 (see below). Using the small-it expansions of £(u) and VKC) 
above we obtain recursively 



An = 2k 2 



A 



1 35 1 25025 

A 3 = iff £;V + 



^2 



S 

Si 

s 2 
#3 



,J_ k 2 u 3 - 667 ^ 2 "< 5 
35 



, )Vr , 5 /,-// ' + 0{ll' , 
Yd k 2 u 2 - f§ fcV + 0(„«) 



■iff k 2 u 3 + 0(u 5 ) 



Up to order 0(u 2 ), only the underlined terms contribute, finally giving 



(u) =1 — k 2 u — ^k^v^lnu 



+ 



1 



-fc 4 (Infc + 7B -|) + - + — j + 



408 



10 35/c 4 715/c 8 



-12a 



(22) 



u 2 + 0(u 3 ). (23) 



2.2 The vacuum contribution 

The T — > limit amounts to neglecting the z 4 terms in Eq. On the k,co plane, below the 
light cone, the solution which is finite for z — > oo and normalised to 1 at z = is 



<j> = \ {k 2 - oo 2 )z 2 K 2 (Vk 2 -uj 2 z). 
Continuing across the light cone into u > k this becomes 

\ m{uj 2 - k 2 )z 2 H^(V^ 2 -k 2 z). 



(24) 
(25) 



The retarded Green's function following from these [S] is given below in (38). 



2.3 The region around the light cone, to w k 

The numerical results presented later show that the dominant contribution comes from the 



region around the light cone. To study this region, we take in ([9|) again u = z as a new 
variable and remove the first derivative term by writing (f)(u) = W(u)Wu/ (1 — u 2 ). The 
equation then becomes (whenever u is used, k = k/(2irT)) 



W"(u) 



(1-u 



2\2 



U + 



k 2 1 



A: 2 



+ 



4u 2 (1 - u 2 ) 



W(u). 



(26) 



To find the behavior for large k, we scale out k by using k 2 l 3 u as the variable. This leads to 
the equation 



1 + 2 — -=— )u + 



k 2 



W"{u) = 

Exactly on the light cone this is simply 

W"(u) = 



+ 3 1 4 

P73~ u + i^ 2 + ^W 3 ' 



W(u). 



(27) 



(28) 
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which is integrable in terms of Airy functions or Bessel functions of order 2/3. The solution 
corresponding to waves infalling into the black hole is given by #2/3 an< ^ the final solution, 
returning to (ft, is 



k^uH^ ku 3 ' 2 ) 



:3 2 / 3 r(§)i 



(29) 



where (-1) 1/3 = § (1 + iy/3). As we shall see concretely, the u 2 term here gives directly the 
Green's function on the light cone. Near the light cone, the term (to 2 — k 2 )/k 2 ~ 2(u — k)/k 
is very small and (oj 2 - k 2 )/k 2 / 3 ~ 2k l / 3 (u) - k). Thus one predicts that the outcome will 
depend mainly on the combination k l / 3 (uj — k). We shall confirm this numerically (see Fig.kl). 



2.4 Method of numerical solution 

The computation of the correlator begins by numerically finding solutions <f>h(z,K), K = 
(u>,k) of (N) representing infalling waves, ~ exp(— icut)(l — z)~ 4tJ//4 . These satisfy 4>* h (z,K) = 
<j>h(z, —K). For u = the situation is particularly simple, one can choose a solution normalised 
to 1, the other is logarithmically divergent. 

Because of the 1/(1 — z 4 ) factor, the integration cannot be started exactly at z = 1. One 
then expands the solution around z = 1. The expansion starts 



<i>h(z) 



(l-z) 



1 + 



(1 - z) (4k 2 - 3w(w - 2i)) 



16 - 8ico 

(1 - z) 2 (16/c 4 - 24A; 2 w 2 + u (9w 3 + 2iu? + 4&j + 32i)) 
128 (oj 2 + 6iu - 8) 



+ ... 



(30) 



up to 20 terms were used. One then expands this solution in terms of the two independent 
solutions at the boundary z = 0, the unnormalisable solution (p u (z, K) and the normalisable 
solution (j> n (z, K): 

<j, h (z, K) = A(K)<j, u (z, K) + B(K)$ n (z, K). (31) 
The expansions start (up to order z 40 were used) as 



(z,K) = 1 + -z 2 (oj 2 - k 2 ) + — z 6 (k 6 



3kW + 3k z (u 4 



>)+0(2 



A: 



2\2 



16 



log(z)cf) n (z,K), 



(32) 



4>n(z,K) 



1 + -Z 2 (k 2 
12 v 



! ) + 



384 



(fc 4 - 2A; 2 a; 2 + w 4 + 192) 



+- 



(A; 6 - 3A; 4 a; 2 + 3A; 2 w 4 + 1344A; 2 



1728w 2 ) 



23040 



+ 0(z 



12 \ 



(33) 



(34) 
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Their Wronskian is 

z 4 ) and W(4> h ,(f> u ) 



t>uK-K4>u = 4z 3 /(l-z 4 ). Similarly, from (|3lJ), W((j) h , 
= -4B{K)z 3 /(l- z A ). 



4A(K)z 3 /(l- 



For the Green's function one needs [6], expanding near z = 0, 
1 <p h (z,K) _u? -k 2 (u 2 



z 3 <t> h {z,K) 



2z 2 



-^(4108 * + 3) + 4 ■ |^ + 0(z 2 ). (35) 



The first two real terms are neglected as contact terms (and they anyway vanish on the light 
cone) and the result for the Green's function is, including 1/(167^5) from the gravity action 



where 



G(K) 



B{K) 



C 3 B{K) 
4vrG 5 A(K) ' 



(36) 



(37) 



Since both Wronskians are ~ z 3 / (1 — z 4 ), the ratio is independent of z and could, in principle, 
be evaluated at any z. In practice, the z-independence is limited by how many terms are 
included in the small-z expansions of 4> u , (fi n . Note that the dimensions of B, (j) n are +4, —4; 
A, 4> u are dimensionless. The neglect of the divergent real terms in the expansion (35) may 
seem somewhat surprising, without counter terms the result (36) only obviously holds for 
the imaginary part. We find that it also reproduces the real part correctly when it can be 
analytically computed. 

As a first application [6] , from the properly analytically continued vacuum solutions ( 24 ) 

(38) 



and (25) one finds that 

G vac (u;, k) 



32 



k 2 ) 2 r -(oj 2 - k 2 ) 

log ~^— A + 2 1E 



Here and until further notice we omit the normalisation factor C 3 / (AttGs) = N 2 /(27r 2 



3 Numerical results for G{k) 

3.1 The static case cu = 



Comparing Eqs. (23) and (31) and changing from k = k/(2nT) back to k = k/(nT) one sees 
that A(k) = 1 and 



B(k) 



i^ 4 [\og\k + lE -\] +i +0(k-*). 



(39) 



There thus is a diverging vacuum contribution plus a constant term 1/10. Including this in 
GVac we write 



G(k) 
G vac (k) 



G V ac(fc) + [G(k) - G vac (k)] = G vac (k) + G T (k), 



± k 4 [log \ k + 7b 



3] , J_ 
4 J ^ 10 



(40) 
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Figure 1: Left panel: Numerically computed finite T part Gr{k) of the correlator, the vacuum part 
(401 is subtracted. Right panel: Check of the analytic expressions (41 1 for the \aige-k behavior of 
Gt(/c). The constants here were originally determined numerically, only later confirmed analytically. 



Furthermore, the large-A: behavior of Gt is 



G T (k) 



48 104448 
35F + 7l5^ + 0( 



12 x 



(41) 



In the numerical integration of ( 13 ) the initial condition is simply <^(1, 0, k) = 1, the other 
solution diverges logarithmically. Because of the 1/(1 — z 4 ) factor the integration is started 
from z = 1 — e' (e ; = 0.2), correcting the initial conditions of <ph an d <ph by the expansion 
(30) (up to (1 — z) 20 was used). Eq. (37) is then evaluated at some z = e (e = 0.2) using the 



expansions (32) and (33) (up to z was used). The finite T part obtained after subtracting 



the vacuum part in ( 40 ) is plotted in Fig. [T] The figure also shows how well its large k behavior 
in Eq. ( |14[ ) is reproduced. One sees that this expansion converges rapidly, already at k = 6 
the error is < 1%. What is important is that there are no terms decreasing more slowly. 

Fig. [T] gives the quantities in dimensionless units. If Gt and k are in physical units, Fig. [T] 
plots G T (fc)/(7rT) 4 vs k/(irT). 



3.2 The case with u/0 

The static correlator G(r) is obtained by Fourier transforming G(k) in the previous subsec- 
tion, but it may be illuminating to see how the same result is obtained with full u dependence, 
using Eq. (§. 



To begin with, from (38) one finds that (u > 0) 

7T 



ImG v 



Pvac(^) 



32 



{u 1 - k 2 ye(oo-k). 



(42) 



Fig.[2]shows how well the numerics approaches this vacuum spectral function. This divergence 
will be subtracted in what follows. 

After vacuum subtraction one finds that p — p vac is non-zero essentially close to the light 



cone. The solution for </> exactly on the light cone is given by the Hankel function in (29) 
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Figure 3: Data for [p(tJ, k) — p vac (uj,k)]/uj, p vac = tt/32 ■ (lu 2 ~ fc 2 ) 2 6(w — k), the left panel plots 



fixed values of k as a function of u>, the right panel the opposite. The red curve shows the result (431 
exactly on the light cone, p(k,k)/k. The red dotted curve is an exponential fit 0.51e 2 ' 5< -' c ~ 10 '. 



The Green's function is evaluated from (37) or directly by computing 2/u ■ <j>' (u) / cft(0) with 
the result 

G(k, k) = (1 + ^) 8 . 6 V3r(2/ 3 ) fc4/3 + w (°- 136092 + i 0-235718) k 4 l\ (43) 



In transforming (29) one must remember that there k = k/(2ivT) while in the numerics k = 
k/(irT). Also note that T(| )T(| ) = 2vr/^/3. Fig. || plots the vacuum-subtracted imaginary 
part either as a function of u for various values of A; or as a function of k for various values 



of u) and one sees that the analytic form ( 43 ) works very well at w = k (note that this is not 
the top of the curves), even down to k>2. The left panel of Fig. [3] shows how the decrease off 
the peak is roughly exponential |6|. For the real part (not shown) the agreement is similar, 
though the A: 4//3 behavior sets in at somewhat larger values of k>10. 

To study in more detail the region around the light cone, Fig. [4] shows the variation when 
one crosses the light cone perpendicularly to it, in the direction of w_ = (jjj — k)/V2 at fixed 
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Figure 4: Data for [p(u>, k) — p vac (uj, k)]/oj when crossing the light cone in a perpendicular direction, 
i.e., in the direction of ui± = (uj ± fc)/\/2, plotted for increasing values of uj + . Note how the peak 
gets narrower with increasing oj + , the curves are of the form fc 1 / 3 ^(A; 1 / 3 ^;-). The real part has a 
similar structure, but the positive peak is at w_ < and the dips at cj_ > are deeper. 




Figure 5: Data for the integral (44) as a function of k. The limits for k — > and k — > oo are 37r/20 



and 47r/20.The right panel has the numerically obtained asymptotic value 0.6286 subtracted. 



ne sees that the prediction made after Eq. u29« is confirmed: the quantity 
in Fig. [4] depends essentially on the combination A; 1 / 3 ^ — k), it is of the form k 1 ^ 3 F(k 1 ^ 3 uj-), 
where -F(O) can be read from (43). Thus while the peak height increases, its width decreases. 



A consequence of this is that the integral over the peak, 



p v£LC (uj,k) 



(44) 



approaches a constant at large k, see Fig. [5| The accuracy is such that within the range 
10. ..80 the last digit in 0.6286 varies between 3. ..8. Subtracting this constant one obtains the 
right panel of Fig. [5] After multiplication by 2/tt this curve is, within numerical accuracy, 
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the same as in the left panel of Fig. [TJ One has, along a rather roundabout way on the k,co 
plane, shown that Eqs. Q and ^ lead to the same result. 

That this numerically obtained constant is actually 7r/5 can be derived as follows by first 
deriving its value at k = 0. A properly subtracted dispersion relation is (at any A:) 

G(w) - G va c(w) - \ G{uj -> oo) - G vac (w -> oo) = / — . (45) 

J-oo vr u/ - u 

For = one can work out the u) — > oo limit of G(uj, k = 0) exactly as the /c — > oo limit of 
G(cj = 0, fc) was worked out in subsection |2.1| with the result G(oo — > oo) — G vac (u; — > oo) 



-3/10 + C(l/w 4 ) [9J. This constant -3/10 is the analogue of the constant 1/10 in Eq. (14) 



If one now takes oj = in (45), the LHS leaves just the constant 3/10 and the RHS is 2/irx 



(44) at k = 0: the value of (44) at k = is 37r/20. Since the k dependence has to be the same 



as that in the static case, Fig. [TJ the value of p4| ) at k = oo is 47r/20. 

4 G(r) 

4.1 The static case u = 

To finally get the correlator in coordinate space we have to compute the integral (the nor- 
malisation £ 3 /(47rG5) will be appended later) 

G(r) = J [GvacW + G T {k)] . (46) 



In the Fourier transform of G va c all terms but the log A; term in (40) are contact terms 



proportional to <5(x) or derivatives thereof, and can be neglected. Using 



\ogk = ^ ^( e -*- e -*'j (47) 

the logfc term becomes 

f t— ^e lk ' x G vac (k) = -oo contact term H \- f — / dk k 5 sm(rk)e~ kt = 

J (2-kY 32fT 2 r J t J 8nr' 

If we introduce /{AitG^) = N^/(27r 2 ), as is appropriate for M = 4 conformal field theory, 
and restore physical units by appropriate powers of Zh, the final result can be written as 

G(r) 15 1 



1 f°° 

N^TY ~ 16*3** + 4^ J dkksHrk)G T (k), f=«Tr. (49) 

Comparing with the next-to-leading order computation in SU(A r c ) Yang-Mills theory [1] 
one notes that the leading UV terms^ are the same if one corrects for an obvious factor of 

lr The leading diagram of the F correlator is reduced to scalar master integrals in Eq. (3.1) of |3j. The 
master integrals are evaluated in Eq. (A. 16) and the final Fourier transformation from G(p) to G(r) is carried 
out using Eq. (5.1) 
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r = 3 T r 



MS 



T= n 



G T (r) 



7rTr 



= 3T, 



-0.0003 - 




NAnT) 



Figure 6: Comparison of the result Gt for the finite T part in ( 49 ) (multiplied by 16 to get a correlator 
of F 2 ) with the NLO computation in SV(N C ) Yang-Mills theory in g] for T = 3T C and T = 12T C 
(dashed). The right panel shows the large k region. The leading 1/r 7 divergence at small r (not shown) 
is the same for the two computations. 




Figure 7: Comparison of the relative magnitudes of the vacuum and finite T parts of (49) (multiplied 
by 16 to get a correlator of F 2 ). Both terms are multiplied by r 7 



" 7 r = nTr. 



16. This arises since [3] studies correlators of F® U F® U while in AdS/CFT duality the scalar 
couples to the operator \ F 2 with the factor -r . 



The finite T part of the full result in (49), multiplied by 16, is compared with the NLO 
computation in [2] in Fig. [6j the leading 1 /r 7 terms are the same. One observes that in the 
finite T part there is a similar region of negative correlator around trTr = 1, but there is one 
crucial difference, the QCD result contains terms ~ (e + p)/r 3 and (e — 3p)/r 3 while there are 
no terms ~ 1/r 3 nor other terms diverging for r — > in the finite T AdS result. The NLO 
QCD result turns positive at irTr ~ 2.6, we do not see this in the AdS result. 

The relative magnitudes of the vacuum and finite T parts of (49) are compared in Fig. [7| 
One sees that the vacuum part dominates for irTr < 1, the finite T part grows in relative 
importance for 1 < nTr < 4 and for nTr > 4 both terms essentially cancel each other. 
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4.2 The case with u^O 

In this case one can write, from Eq. ([5]), 



G{r) 



1 



2n 2 r 
1 



dk k sin(rfe) 



7T 



2n 2 r 7r 7 



1 f 15vr 2 



dfc sin(rfc) x 

rD0 7T 



(a; 2 



+ 



a; 32 
d/c fe sin(r/j) 



du p{uj,k) - p vauC (uj,k) 



7T 







+ 

g - Pvac 



Pvac _ Q.6286 + q >6286 



CO 



0.6286 



+ const • <5(x). (50) 



Here the Fourier transform of the vacuum contribution is computed by introducing a conver- 
gence improving factor e~™ . The vacuum integrals over u and k can then be done analytically 
and the limit t — > leads to the expression in (50) above. Since the square bracket in the 
integrand in the RHS, after multiplication with 2/tt, coincides with Gx(k), this result is the 
same as in (49). 



5 G(r,k = 0) and G(r = 0,r) 



Various r dependent correlators can also be computed using the spectral representation Q. 
A particularly simple case is the t dependent correlator at zero spatial momentum, obtained 
by summing over spatial volume. Using the spectral representation Q and the explicit form 
(42) one has, within the range < r < 1/(2T), 



G(r,k = 0) 



duj 



: : !2 



w 4 + i( / 9(w,0)-p vac (w,0)) 



cosh[(l-2Tr)| Cj] 
sinh 2 0J 



4vr 5 



[((5,Tt) + ((5,1 -Tt)] + G 2 (Tt) 



(51) 



where one has noted that in the dimensionless units used here (3oj — > ituj and where G2 has to 
be integrated numerically. The integrand is similar to the k = 0.5 curve in Fig. [3j multiplied 
by u. For small x, 



C(5,x) + C(5,l 



1 



+ 2C(5) + 30C(7)x 2 + O(x 4 



(52) 



Physical dimensions of G are restored by multiplying by (ttT) 5 so that the first term produces 
the T independent UV divergence 3/(4r 5 ). Again, there are no further divergent terms. One 
finds that the component G2 is numerically insignificant even at TV = ^ ; where it is largest 
relative to the £ function terms. 

A more complicated case is G(r,r). We give, for completeness, only the part obtained by 
inserting to @ the vacuum spectral function d42|): 



G V s, c (t, r) 



7r 2 r ! 



E 



+ 



7r 2 (n + Tr) 2 /r 2 + l] 4 [vr 2 (n + 1 - Tr) 2 /r 2 + l] 4 



(53) 
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The sum can be expressed in terms of derivatives of the gamma function, but this is not very 
illuminating. For equal r correlators one has 



G vac (r = 0,r) 



6 



— r 4 csch 4 (r) + — r 4 coth 2 (r)csch 2 (r) + -r 3 coth(r)csch 2 (r) + 



— r 2 csch 2 (r) H rcoth(r) 

16 W 16 W 



(54) 



where the quantity in the brackets = 1 + r 8 /4725 + 0(r 10 ). Including the factor N 2 /{2ir 2 ) 
and multiplying by 16 to get correlators of F 2 , the leading singular term ~ 48iV 2 /(7r 4 r 8 ) is 
seen to agree with that in [3]. 



6 Conclusions 

Motivated by needs of lattice work [1], we have in this article computed the r dependence 
of the r integrated finite temperature correlator (F 2 (t, x)i ?2 (0, 0)) and compared it with 
the same in SU(A r c ) Yang-Mills theory [3J. Our computation based on AdS/CFT duality 
applies to N = 4 conformal strongly coupled theory and there is no reason for the results to 
coincide. The leading UV vacuum part ~ 1/t 7 , t — > 0, coincides since it is independent of 
the coupling. However, the leading finite T parts differ, in the QCD-like case there are terms 
~ (e + p)/r 3 and ~ (e — 3p)/r 3 , while for AdS/CFT we find no terms ~ T 4 /r 3 nor other 
finite T terms diverging for r — > 0. Gt is just some function ^^(vrTr), finite for r — > 0. In 
the range r ~ few l/(irT) both finite T correlators are negative. For r > 4/(ttT)) the AdS 
finite T correlator essentially cancels the vacuum one. 

The computation was carried out both taking u) = from the outset, but also by integrating 
the weighted spectral function over uj. These results coincide due to observed structure near 
the light cone, ui = k. It would be interesting to have analytic control of observed scaling as 
a function of k l /^(uj — k). 

An obvious task for the future would be carrying out a similar computation for holographic 
QCD models [T2] breaking conformal invariance. This will also shed light on the difference 
between (F 2 (x)F 2 (0)) and {F 2 (x)F 2 (0)) correlators. 
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